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The degree to which a pure quantum state is entangled can be characterized by the distance or 
angle to the nearest unentangled state. This geometric measure of entanglement, already present 
in a number of settings (see Shimony [[ and Barnum and Linden Q), is explored for bipartite 
and multipartite pure and mixed states. The measure is determined analytically for arbitrary two- 
qubit mixed states and for generalized Werner and isotropic states, and is also applied to certain 
multipartite mixed states. In particular, a detailed analysis is given for arbitrary mixtures of three- 
qubit GHZ, W and inverted-W states. Along the way, we point out connections of the geometric 
measure of entanglement with entanglement witnesses and with the Hartree approximation method. 
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I. INTRODUCTION 

Only recently, after more than half a century of exis- 
tence, has the notion of entanglement become recognized 
as central to quantum information processing |3(. As 
a result, the task of characterizing and quantifying en- 
tanglement has emerged as one of the prominent themes 
of quantum information theory. There have been many 
achievements in this direction, primarily in the setting 
of bipartite systems 0]. Among these, one highlight is 
Wootters' formula [5| for the entanglement of forma- 
tion for two-qubit mixed states; others include corre- 
sponding results for highly symmetrical states of higher- 
dimensional systems 0, 0|- The issue of entanglement 
for multipartite states poses an even greater challenge, 
and there have been correspondingly fewer achievements: 
notable examples include applications of the relative en- 
tropy p|, negativity and Schmidt measure [Icj . 

In this Paper, we present an attempt to face this chal- 
lenge by developing and investigating a certain geomet- 
ric measure of entanglement (GME), first introduced by 
Shimony in the setting of bipartite pure states and 
generalized to the multipartite setting (via projection 
operators of various ranks) by Barnum and Linden 0. 
We begin by examining this geometric measure in pure- 
state settings and establishing a connection with entan- 
glement witnesses, and then extend the measure to mixed 
states, showing that it satisfies certain criteria required 
of good entanglement measures. We demonstrate that 
this geometric measure is no harder to compute than 
the entanglement of formation E-p, and exemplify this 
fact by giving formulas corresponding to Ep for (i) ar- 
bitrary two-qubit mixed, (ii) generalized Werner, and 
(iii) isotropic states. We conclude by applying the geo- 
metric entanglement measure to certain families of mul- 
tipartite mixed states, for which we provide a practi- 
cal method for computing entanglement, and illustrate 
this method via several examples. In particular, a de- 
tailed analysis is given for arbitrary mixture of three- 
qubit GHZ, W, and inverted-W states. 



It is not our intention to cast aspersions on exisiting 
approaches to entanglement; rather we simply wish to 
add one further element to the discussion. Our discus- 
sion focuses on quantifying multipartite entanglement in 
terms of a single number rather than characterizing it. 

The structure of the paper is as follows. In Sec. [H] 
we describe the basic geometric ideas on quantifying en- 
tanglement geometrically in the setting of pure quan- 
tum states, and establish a connection with the Hartree 
approximation method and entanglement witnesses. In 
Sec. IIIII we extend the definition of the GME to mixed 
states, and show that it is an entanglement monotone. 
In Sec. IIVI we examine the GME for several families of 
mixed states of bipartite systems: (i) arbitrary two-qubit 
mixed, (ii) generalized Werner and (iii) isotropic states 
in bipartite systems, as well as (iv) certain mixtures of 
multipartite symmetric states. In Sec. ^ we give a de- 
tailed application of the GME to arbitrary mixtures of 
three-qubit GHZ, W and inverted-W states. In Sec. IVII 
we discuss some open questions and further directions. 
In App. |l|we briefly review the Vollbrecht- Werner tech- 
nique used in the present Paper. 



II. BASIC GEOMETRIC IDEAS AND 
APPLICATION TO PURE STATES 

We begin with an examination of entangled pure 
states, and of how one might quantify their entanglement 
by making use of simple ideas of Hilbert space geometry. 
Let us start by developing a quite general formulation, 
appropriate for multipartite systems comprising n parts, 
in which each part can have a distinct Hilbert space. 
Consider a general n-partite pure state 

W= E Xp lP2 ---pJe$e$---e<£). (1) 

Pl—Pn 
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One can envisage a geometric definition of its entangle- 
ment content via the distance 

d = min|||V.'}- |0)|| (2) 

m 

between \ip) and the nearest separable state |0) (or equiv- 
alently the angle between them). Here, |0) = \(jr- 1 ') 
is an arbitrary separable (i.e., Hartree) n-partite pure 
state, the index i = 1 . . . n labels the parts, and 

|^))=£c«|e«). (3) 

Pi 

It seems natural to assert that the more entangled a state 
is, the further away it will be from its best unentangled 
approximant (and the wider will be the angle between 
them). 

To actually find the nearest separable state, it is con- 
venient to minimize, instead of d, the quantity 

lll^-l^ll 2 , (4) 

subject to the constraint ((f>\<j>) = 1. In fact, in solv- 
ing the resulting stationarity condition one may restrict 
one's attention to the subset of solutions \<j>) that obey 
the further condition that for each factor \4>^) one has 
(0W|0«) = 1. Thus, one arrives at the nonlinear eigen- 
problem for the stationary |0): 

V v* ••■r (i) ...>) - A f W' 

Xpip2---pn C J>i ' ' ' C Pi ' ' ' C J)j — ^ C J}} ! (Pb) 

Pl---Pi---Pn 

where the eigenvalue A is associated with the Lagrange 
multiplier enforcing the constraint (0|0) — 1, and the 
symbol " denotes exclusion. In basis-independent form, 
Eqs. © read 

® \4>W)) = A(0«|, (6a) 

( ® (0 (J) |W> = A|0«). (6b) 

From Eqs. JSJ) or © one readily sees that the eigenvalues 
A are real, in [—1,1], and independent of the choice of 

the local basis {|e p ^)}. Hence, the spectrum A is the 
cosine of the angle between and |0); the largest, A max , 
which we call the entanglement eigenvalue, corresponds 
to the closest separable state and is equal to the maximal 
overlap 

A max = max 11(01-0)11, (7) 

4> 

where \<p) is an arbitrary separable pure state. 

Although, in determining the closest separable state, 
we have used the squared distance between the states, 
there are alternative (basis-independent) candidates for 



entanglement measures which are related to it in an ele- 
mentary way: the distance, the sine, or the sine squared 
of the angle 9 between them (with cos 6 = Re (ip\<j))). We 
shall adopt E sin 2 = 1— A^ ax as our entanglement measure 
because, as we shall see, when generalizing E sin 2 to mixed 
states we have been able to show that it satisfies a set of 
criteria demanded of entanglement measures. We remark 
that determining the entanglement of \ip) is equivalent to 
finding the Hartree approximation to the ground-state of 
the auxiliary Hamiltonian Ti = — \ijj){ip\ [Tj. 

In bipartite applications, the eigenproblem (JSJ is in 
fact linear, and solving it is actually equivalent to find- 
ing the Schmidt decomposition 0. Moreover, the en- 
tanglement eigenvalue is equal to the maximal Schmidt 
coefficient. By constrast, for the case of three or more 
parts, the eigenproblem is a nonlinear one. As such, one 
can in general only address it directly, i.e., by determin- 
ing the eigenvalues and eigenvectors simultaneously and 
numerically. Yet, as we shall illustrate shortly, there do 
exist certain families of pure states whose entanglement 
eigenvalues can be determined analytically. 

A. Illustrative examples 

Suppose we are already in possession of the Schmidt 
decompostion of some two-qubit pure state: 

IV;} = Vp|oo) + v / i _ p|ii>. (8) 

Then we can read off the entanglement eigenvalue: 

A max = maxjy^, V 1 "Pi- ( 9 ) 

Now, recall that the concurrence C for this state is 
2y/p(l — p). Hence, one has 

AL X = \ (i + Vi - c 2 ) , (io) 

which holds for arbitrary two-qubit pure states. 

The possession of symmetry by a state can alleviate 
the difficulty associated with solving the nonlinear eigen- 
problem. To see this, consider a state 

W)= E X Pl p 2 ---pJe£e$---eW) (11) 

Pl—Pn 

that obeys the symmetry that the nonzero amplitudes x 
are invariant under permutations. What we mean by this 
is that, regardless of the dimensions of the factor Hilbert 
spaces, the amplitudes are only nonzero when the in- 
dices take on the first v values (or can be arranged to do 
so by appropriate relabeling of the basis in each factor) 
and, moreover, that these amplitudes are invariant under 
permutations of the parties, i.e., Xa-xa- 2 ---cr„ = XpxP2--Pni 
where the cr's are any permutation of the p's. (This 
symmetry may be obscured by arbitrary local unitary 
transformations.) For such states, it seems reasonable 
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to anticipate that the closest Hartree approximant re- 
tains this permutation symmetry. Assuming this to be 
the case — and numerical experiments of ours support this 
assumption — in the task of determining the entanglement 
eigenvalue one can start with the Ansatz that the closest 
separable state has the form 

10) fe c il e f)j . ( 12 ) 

i.e., is expressed in terms of copies of a single factor state, 

(i) 

for which c - = Cj. To obtain the entanglement eigen- 
value it is thus only necessary to maximize Re (cj)\ip) with 
respect to {cj}j =1 , a simpler task than maximization over 
the y™—! di amplitudes of a generic product state. 

To illustrate this symmetry-induced simplification, 
we consider several examples involving permutation- 
invariant states, first restricting our attention to the case 
v = 2. The most natural realizations are n-qubit sys- 
tems. One can label these symmetric states according to 
the number of O's, as follows [l2^ : 

| 5 (n,fc)) = y5EM £ I0__01__l). (13) 

permutations y. n—k 

As the amplitudes are all positive, one can assume that 
the closest Hartree state is of the form 

\4>) = (^>n + ^T—p\\)f\ (w) 

for which the maximal overlap (w.r.t. p) gives the entan- 
glement eigenvalue for |S(n, k))\ 

W»,*) = {^^{^ (^) 2 • W 

For fixed n, the minimum A max (and hence the maximum 
entanglement) among the |S(n, fc))'s occurs for k = n/2 
(for n even) and k = (n ± l)/2 (for n odd). In fact, for 
fixed n the general permutation-invariant state can be ex- 
pressed as J2 k ak\S(n, k)) with J2 k \ctk\ 2 — 1- The entan- 
glement of such states can be addressed via the strategy 
that we have been discussing, i.e., via the maximization 
of a function of (at most) three real parameters. The 
simplest example is provided by the 71GHZ state: 

|nGHZ) = (|S(n, 0)) + |S(n, n))) /V2. (16) 

It is easy to show that (for all n) A max (nGHZ) = 1/V2 
and E sin i = 1/2. 

We now focus our attention on three-qubit settings. Of 
these, the states |5(3,0)) = |000) and 15(3,3)) = |111) 
are not entangled and are, respectively, the components 
of the the 3-GHZ state: |GHZ) = (|000) + |111>)/V2. 
The states |S(3,2)) and |S(3,1)), denoted 

|W) = |S(3,2)) = (|001) + |010) + |100))/V3, (17a) 
|W) = |S(3,1)) = (|110) + |101) + |011))/V3, (17b) 
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FIG. 1: Entanglement of the pure state a/s | W) + y/l - s | W) 
vs. s. This also turns out to be the entanglement curve for 
the mixed state s |W)(W| + (1 - s)|W)(W|. 

are equally entangled, having A max = 2/3 and -E S i n 2 = 
5/9. 

Next, consider a superposition of the W and W states: 

|WW(s») = y^|W) + vT^ie^lW). (18) 

It is easy to see that its entanglement is independent of 
<j>: the transformation (|0), |1)} -> (|0),e _i *|l)} induces 

|WW(s, 0)) e"^|WW(s,0)). To calculate A max , we 
assume that the separable state is (cos#|0) + sin^jl))® 3 
and maximize its overlap with |WW(s, 0)). Thus we find 
that the tangent t = tan 9 is the particular root of the 
polynomial equation 

VI -st 3 + 2Vst 2 - 2\/l - st - Vs = (19) 

that lies in the range t G [a/1/2, y/2\. Via 9(s), A max 
(and thus E sin 2 = 1 — A max ) can be expressed as 

A max (s) = - [Vs cos 6(s) + Vl - s sin 0(a)] sin26»(s). (20) 

In Fig.^ w e show E sin 2 (|WW(s, </>))) vs. s. In fact, A max 
for the more general superposition 

|SS„ ;felfc2 (r, 4>)) = |S(n, h)) + VT^e^ \S(n, fc 2 )) 

(21) 

(with k\ 7^ k%) turns out to be independent of 0, as in 
the case of | WW(s, (/>)), and can be computed in the same 
way. We note that although the curve in Fig.Q]is convex, 
convexity does not hold uniformly over fci and hi. 

As our last pure-state example in the qubit setting, we 
consider superpositions of W and GHZ states: 

|GW(s,0)) = V^|GHZ) + Vl — ^e^|W). (22) 

For these, the phase <j> cannot be "gauged" away and, 
hence, E sin 2 depends on <f>. In Fig. [5] we show E sin 2 vs. s 
at 4> = and <f> — ir (i.e., the bounding curves), as well 
as -Esin 2 for randomly generated values of s G [0, 1] and 
4> S [0, 2tt] (dots). It is interesting to observe that the 
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FIG. 2: Entanglement of |GW(s, </))) vs. s. The upper curve is 
for <f> = 7T whereas the lower one is for (f> = 0. Dots represent 
states with randomly generated s and </>. 



V state has higher entanglement than the '0' does. As 
the numerical results suggest, the (0-parametrized) E sin 2 
vs. s curves of the states \GW(s,(j>)} lie between the c 7r' 
and '0' curves. 

We remark that, more generally, systems comprising n 
parts, each a d-level system, the symmetric state 



feo fci 



with ki = n, has entanglement eigenvalue 

d-l 



A max (n; {k}) 



n 

i=0 



(23) 



(24) 



One can also consider other symmetries. For example, 
for the totally antisymmetric (viz. determinant) state of 
n parts, each with n levels, 



|Dct„> 



n 



«1, 



(25) 



i\ ,...,i n =l 



it has been shown by Bravyi |13l | that the maximal 
squared overlap is A^ ax = l/n\. Bravyi also generalized 
the anti-symmetric state to the n = p cP-partite determi- 
nant state via 



0(2) 



(0,0,. ..,0,0), 
(0,0,. ..,0,1), 



(j){dP-l) = (d-l,d-l,...,d-l,d-2), 
') = (d-l, d-l,..., d-l, d-l), 



and 



l p etn,d) = j== e ilt ... ! i dP \(f)(ix),...,(j>(i (t p)). 



In this case, one can show that A 2 lax = [(e? p )!] 1 



(26) 



B. Connection with entanglement witnesses 

We now digress to discuss the relationship between the 
geometric measure of entanglement and another entan- 
glement property — entanglement witnesses. The entan- 
glement witness W for an entangled state p is defined to 
be an operator that is (a) Hcrmitian and (b) obeys the 
following conditions 14]: 

(i) Tr(yVer) > for all separable states a, and 

(ii) Tr(Wp) < 0. 

Here, we wish to establish a correspondence between 
A max for the entangled pure state and the optimal 
element of the set of entanglement witnesses W for 
that have the specific form 



W 



A 2 1 



mm, 



(27) 



this set being parametrized by the real, non-negative 
number A 2 . By optimal we mean that, for this spe- 
cific form of witnesses, the value of the "detector" 
Tr(W\ip}(ip\) is as negative as can be. 

In order to satisfy condition (i) we must ensure that, 
for any separable state a, we have Tr(Wo-) > 0. As 
the density matrix for any separable state can be de- 
composed into a mixture of separable pure states [i.e., 
(j = J^i where {|</>i)} are separable pure states], 

condition (i) will be satisfied as long as Tr(W|</>)(0|) > 
for all separable pure states \<j)). This condition is equiv- 
alent to 

- IIW#)H 2 > (for all separable \<j>)), (28) 
which leads to 

A 2 >max||(^|0)|| 2 =A 2 nax (|^)). (29) 

m 

Condition (ii) requires that Ti (W\ip)(ip\) < 0, in order 
for KV to be a valid entanglement witness for this 
gives A 2 — 1 < 0. Thus, we have established the range of 
A for which A 2 1 — is a valid entanglement witness 

for lib): 



< 1. 



(30) 



With these preliminaries in place, we can now estab- 
lish the connection we are seeking. Of the specific fam- 
ily (|27|l of entanglement witnesses for \ip), the one of the 
form Wopt = A 2 nax (|?/'))1 — ls optimal, in the sense 

that it achieves the most negative value for the detector 
Tr(W pt|V)(^l): 

minTr(nWM) = Tr(W op t|^)^|) = 

(31) 

where W runs over the class <|27|) of witnesses. 

We now look at some examples. For the GHZ state 
the optimal witness is 



Wghz = ^1- | GHZ) (GHZ | 



(32) 
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and Tr(W G Hz|GHZ)(GHZ|) = -E sin 2 (|GHZ)) = -1/2. 
Similarly, for the W and inverted-W states we have 

W W = ^1-|W)(W| and W^ = ^1-|W)(W| (33) 

and Tr(W w |W)(W|) = -E sin 2(\W)) = -5/9, and simi- 
larly for |W). For the four-qubit state 

|*) = (|0011) + |0101) + |0110) + |1001) + |1010) + |1100))/V6 

(34) 

the optimal witness is 

W* = fl-|*)(*l (35) 
o 

and Ti(W*\*)m) = -£w(|*)) - -5/8. 

Although the observations we have made in this sec- 
tion are, from a technical standpoint, elementary, we nev- 
ertheless find it intriguing that two distinct aspects of 
entanglement — the geometric measure of entanglement 
and entanglement witnesses — are so closely related. Fur- 
thermore, this connection sheds new light on the content 
of the geometric measure of entanglement. In particular, 
as entanglement witnesses are Hermitian operators, they 
can, at least in principle, be realized and measured lo- 
cally [3. Their connection with the geometric measure 
of entanglement ensures that the geometric measure of 
entanglement can, at least in principle, be verified exper- 
imentally. 

III. EXTENSION TO MIXED STATES 

The extension of the GME to mixed states p can be 
made via the use of the convex roof (or hull) construction 
[indicated by "co"], as was done for the entanglement of 
formation (see, e.g., Ref. 0]). The essence is a minimiza- 
tion over all decompositions p = Y^iPi m to pure 
states, i.e., 

E{p) = (coE puIe )(p) = min V" p t -Epure(IV'i))- (36) 
{p<>iM — 1 

Now, any good entanglement measure E should, at least, 
satisfy the following criteria (c.f. Refs. [Si llfiL Il7j ) : 

CI. (a) E{p)>0- (b) E(p)=0 if p is not entangled. 

C2. Local unitary transformations do not change E. 

C3. Local operations and classical communication 
(LOCC) (as well as post-selection) do not increase 
the expectation value of E. 

C4. Entanglement is convex under the discarding of in- 
formation, i.e., J2iPi E (pi) > E(J2iPiPi)- 

The issue of the desirability of additional features, such as 
continuity and additivity, requires further investigation, 
but C1-C4 are regarded as the minimal set, if one is to 
guarantee that one has an entanglement monotone 0. 



Does the geometric measure of entanglement obey Cl- 
4? From the definition 136( 1 it is evident that CI and 
C2 are satisfied, provided that E puic satisfies them, as it 
does for -E pU re being any function of A max consistent with 
CI. It is straightforward to check that C4 holds, by the 
convex hull construction. The consideration of C3 seems 
to be more delicate. The reason is that our analysis of 
whether or not it holds depends on the explicit form of 
E puIe . For C3 to hold, it is sufficient to show that the 
average entanglement is non-increasing under any trace- 
preserving, unilocal operation: p — > Ylk^kpVk' where 

(i) 

the Kraus operator has the form V k = 1 ® • • • 1 <g> V k y <8> 

1 • • • (g) 1 and Ylk ^fe ^fc = ft- Furthermore, it suffices to 
show that C3 holds for the case of a pure initial state, 
i.e., p = |V>)(# 

We now prove that for the particular (and by no means 
un-natural) choice -E pur c = E sin 2 , C3 holds. To be pre- 
cise, for any quantum operation on a pure initial state, 
i.e., 

l^(^hE/*H^)!. ( 37 ) 

we aim to show that 

J> fe £w (Yk\i>)/^Pk)<Ean'(m, (38) 
k 

where p k = Tr Vk\ij) (ip\V^ = (ip\V^V k \ip), regardless of 
whether the operation {14} is state-to-state or state-to- 
ensemble. Let us respectively denote by A and A*, the 
entanglement eigenvalues corresponding to |^») and the 
(normalized) pure state Vk\i>) I \fp~k~ ■ Then our task is to 
show that J2kP k ^1 — °f which the left hand side is, 
by the definition of A k , equivalent to 




. (39) 

Without loss of generality, we may assume that it is the 
first party who performs the operation. Recall that the 
condition © for the closest separable state 

|0) = |fi)i <g> |7>2... n (40) 

can be recast as 

2 ...„(7|^)i...„=A|a) 1 . (41) 
Then, by making the specific choice 

(&l = ((fi|^ (1)t A/ffl0®<7l, (42) 
where q k = (a|V^ V^ 1 ^), we have the sought result 
]>> fe A 2 fe = ]r max ||<&|W>|| 2 

k k 

> A 2 ^((5|y«VW|5)/V*) 2 = A 2 . (43) 
k 
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Hence, the form 1 — A 2 , when generalized to mixed states, 
is an entanglement monotone. We note that a differ- 
ent approach to establishing this result has been used by 
Barnum and Linden [2(. Moreover, using the result that 
^ fe PfeA| > A 2 , one can further show that for any convex 
increasing function f c (x) with x € [0, 1], 



(44) 



Therefore, the quantity const. — f c (A 2 ) (where the const. 
is to ensure the whole expression is non-negative), when 
extended to mixed states, is also an entanglement mono- 
tone, hence a good entanglement measure. For the fol- 
lowing discussion we simply take E = 1 — A 2 . 



IV. ANALYTIC RESULTS FOR MIXED STATES 



Before moving on to the terra incognita of mixed mul- 
tipartite entanglement, we test the geometric approach in 
the setting of mixed bipartite states, by computing E sin 2 
for three classes of states for which Ep is known. 



A. Arbitrary two-qubit mixed states 



For these we show that 



(45) 



where C(p) is the Wootters concurrence of the state p. 
Recall that in his derivation of the formula for Ep , Woot- 
ters showed that there exists an optimal decomposition 
p = J2iPi m which every has the concurrence 

of p itself. (More explicitly, every has the identical 
concurrence, that concurrence being the infimum over all 
decompositions.) By using Eq. (|10|l one can, via Eq. 1)45(1 . 
relate E sin 2 to C for any two-qubit pure state. As E s - ln 2 is 
a monotonically increasing convex function of C € [0, 1], 
the optimal decomposition for -E s j n 2 is identical to that 
for the entanglement of formation Ep. Thus, we see that 
Eq. H45|) holds for any two-qubit mixed state. 

The fact that E s[n 2 is related to Ep via the concur- 
rence C is inevitable for two-qubit systems, as both are 
fully determined by the one independent Schmidt coeffi- 
cient. We note that Vidal has derived this expression 
when he considered the probability of success for con- 
verting a single copy of some pure state into the desired 
mixed state, which gives a physical interpretation of the 
geometric measure of entanglement. Unfortunately, this 
connection only holds for two-qubit states. 



B. Generalized Werner states 

Any state pw of a C d <g> C d system is called a general- 
ized Werner state if it is invariant under 



Pi :p 



dU(U <E)U)p(U^ <E>U^), (46) 



where U is any element of the unitary group U(d) and 
dU is the corresponding normalized Haar measure. Such 
states can be expressed as a linear combination of two 
operators: the identity 1, and the swap F = ■ \ij){ji\, 

i.e., pw = al + 6F, where a and b are real parameters 
related via the constraint Trpw = 1 • This one-parameter 
family of states can be neatly expressed in terms of the 
single parameter / = Tr(pwF): 



Pw(f) = 



d 2 - fd 
d 4 -d 2 



a ® a 



fd 2 - d , 
d 4 -d 2 



(47) 



By applying to -E sin 2 the technique by developed by Voll- 
brecht and Werner for Ep(pw) [see Ref. [6] or App. |A"| . 
one arrives at the geometric entanglement function for 
Werner states: 

1 



^W/D^^-VKF) for/<0, (48) 

and zero otherwise. 

The essential points of the derivation are as follows: 
(i) In order to find the set M pw (see App. it is suffi- 
cient, due to the invariance of pw under Pi, to consider 



any pure state |4>) = J2jk®jk\e 



(i) 



,(2) 



that has a 



diagonal reduced density matrix Tr2 1 *&) (*& | and the value 
Tr(|<f>) (<f>|_F) equal to the parameter / associated with 
the Werner state pw(/)- It can be shown that 



£ sin KI*><*l)>~(l-A/l-(/-E A «) 2 



(49) 



where Xu = \§a\ 2 . 

(ii) If / > 0, we can set the only nonzero elements of 
|$) to be $ i2 , . . ., . . ., $id such that |$ i4 | 2 = /, 
this state obviously being separable. Hence, for / > 
we have E sin 2(p w (f)) = 0. On the other hand, if / < 
then any nonzero A^; would increase (/ — X)i^") 2 and, 
hence, increase the value of E(\Q) not conforming 
with the convex hull. Thus, for a fixed value of /, the 
lowest possible value of the entanglement E(\&) ($|) that 
can be achieved occurs when — and there are only 
two nonzero elements $y and (i ^ j). This leads to 



min E(m($\) 

|$) at fixed / 



|(i-v^7 



Thus, as a function of /, e(/) is given by 



\(l-^T~p) for/ 



<0, 



for / > 0, 



(50) 



(51) 



which, being convex for / S [—1, 1], gives the entangle- 
ment function l|48|l for Werner states. 
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FIG. 3: Entanglement curve for the mixed state p7-2,s{r) (full 
line) constructed as the convex hull of the curve for the pure 
state |SS7 ; 2,5(V, </))) (dashed in the middle; full at the edges). 



C. Isotropic states 

These are states invariant under 

P 2 : p-> (dU(U® U*)p(U^ ®U*^), (52) 



and can be expressed as 



1-F 



Piso (F) = -^-^ (i - |$ + >(<f + |) + F|$ + >($ + |, (53) 

where |$+) = ^Ei=iN) and F e t ' 1 ]- For F e 
[0, 1/d] this state is known to be separable [19|]. By fol- 
lowing arguments similar to those applied by Terhal and 
Vollbrecht for i?F(piso) one arrives at 

E sin 2 ( Piso (F)) = 1 - ^(VF+y/Q. -F)(d~l)) 2 , (54) 

for F > 1 /d. The essential point of the derivation is the 
following Lemma (cf. Ref. Q): 

Lemma 1 . The entanglement E sin 2 for isotropic states in 
C d (g) C d for F € 1] is given by 

£w(p iso (iO)=co(E(F)), (55) 
where co(i?(F)) is the convex hull of the function R and 



D. Mixtures of multipartite symmetric states 

Before exploring more general mixed states, it is useful 
to first examine states with high symmetry. With this in 
mind, we consider states formed by mixing two distinct 
symmetric states (i.e., k\ ^ fo): 



Pn-MkJc) = r \ S ( n > k l)) ( S ( n > fc l)l 

+(l-r)|S(n, fc 2 ))(S(n,fc 2 )|. 



(58) 



From the independence of E sin 2 (\SS n ;k 1 k 2 (' r , 4>)}) on 4> 
and the fact that the mixed state p n -k 1 k 2 ( r ) is invariant 
under the projection 



2ir 



(59) 



with U : {|0),|1)} -> {|0>,e-^|l>}, we have that 
^sin 2 {Pn;k 1 k 2 ( r )) vs - r can be constructed from the con- 
vex hull of the entanglement function of ISSn^kj (r, 0)) 
vs. r. An example, (n, fcijfe) = (7,2,5), is shown in 
Fig. |3 If the dependence of -E s i n 2 on r is already con- 
vex for the pure state, its mixed-state counterpart has 
precisely the same dependence. Figure ^ f° r which 
(n, fcijfe) = (3,1,2), exemplifies such behavior. More 
generally, one can consider mixed states of the form 



p({p}) = ^2p k \S(n,k))(S(n,k)\. 



(60) 



The entanglement -E m ixcd({p}) can then be obtained as 
a function of the mixture {p} from the convex hull of 
the entanglement function E pnre ({q}) for the pure state 
J2k Tnat is > E mix({p}) = CoE purc ({q} = 

{p}). Therefore, the entanglement for a mixture of sym- 
metric states \S(n, k)) is known, up to some convexifica- 
tion. 



APPLICATION TO ARBITRARY MIXTURE 
OF GHZ, W AND INVERTED-W STATES 



Having warmed up in Sec. lIVDl bv analyzing mixtures 
of multipartite symmetric states, we now turn our atten- 
tion to mixtures of three-qubit GHZ, W and inverted- W 
states. 



R(F) = l-max{ Ml |F = (^V^l) 2 M ^IM = l}- 

i=l i=l 

(56) 

Straightforward extremization shows that 



R(F) = 1 




(57) 



which is convex, and hence co(R(F)) = R(F). Thus 
we arrive at the entanglement result for isotropic states 
given in Eq. Q54)). 



A. Symmetry and entanglement preliminaries 

These states are important, in the sense that all pure 
states can, under stochastic LOCC, be transformed either 
to GHZ or W (equivalently inverted- W) states. It is thus 
interesting to determine the entanglement content (using 
any measure of entanglement) for mixed states of the 
form: 



p(x, y) = x|GHZ) (GHZ| + y| W) (W| + (1 -x - 



y)|W)(W|, 
(61) 
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FIG. 4: Entanglement vs. the composition of the pure state 
\ip(x,y)). This entanglement surface is not convex near 
(x,y) — (1,0), although not obvious from the plot. 

where x,y > and x + y < 1. This family of mixed 
states is not contained in the family H6()|l . as |GHZ) = 
(|S(3,0)) + |S(3,3)))/-\/2. The property of p(x,y) that 
facilitates the computation of its entanglement is a cer- 
tain invariance, which we now describe. Consider the 
local unitary transformation on a single qubit: 

|0) - |0), (62a) 
11} - 9 fc |l>, (62b) 

where g = exp(27ri/3), i.e., a relative phase shift. This 
transformation, when applied simultaneously to all three 
qubits, is denoted by Uk- It is straightforward to see that 
p{x, y) is invariant under the mapping 

1 3 

T> A :p^-Y,U kP Ul (63) 
6 k=i 

Thus, we can apply Vollbrecht- Werner technique [6j to 
the compution of the entanglement of p(x, y). 

Now, the Vollbrecht- Werner procedure requires one to 
characterize the set S- lm of all pure states that are in- 
variant under the projection P4. Then, the convex hull 
of E s - ln 2 (p) need only be taken over Si nv , instead of the 
set of all pure states. However, as the state p(x,y) is 
a mixture of three orthogonal pure states (|GHZ), |W) 
and |W)) that are themselves invariant under P4, the 
pure states that can enter any possible decomposition of 
p must be of the restricted form: 

a|GHZ) +/3|W) +7|W), (64) 

with |a| 2 + |/3| 2 + |7| 2 = 1. Thus, there is no need to char- 
acterize S'inv, but rather to characterize the pure states 
that, under P4, are projected to p(x, y). These states are 
readily seen to be of the form: 




X 



FIG. 5: Entanglement of the pure state \ip(x,y = (1 — 

x)/2)> = v ^|GHZ) + v'(l-x)/2|W> + V(l-x)/2|W> vs. x. 
This shows the entanglement along the diagonal boundary 
x + 2y = 1. Note the absence of convexity near x = 1; this 
region is repeated in the inset. 

Of these, the least entangled state, for (x, y) given, has 
all coefficients non- negative (up to a global phase), i.e., 

\i/>(x,y)) = ^|GHZ) + Vy|W) + v / l-x-y|W). (66) 

The entanglement eigenvalue of \tp(x, y)) can then be 
readily calculated, and one obtains 

A(^^)= (1+ 1 t2)| {/f( 1+t3 ) + V^t+^/3(l-x-y)t 2 

' ' " ' (67) 
where t is the (unique) non-negative real root of the fol- 
lowing third-order polynomial equation: 

+ v/3(l -x-y){-t 3 + 2t) = 0. (68) 

Hence, the entanglement function for \ip(x,y)}, i.e., 
Etp(x, y) = 1—A(x, y) 2 , is determined (up to the straight- 
forward task of root- finding). 

B. Finding the convex hull 

Recall that our aim is to determine the entanglement 
of the mixed state p(x,y). As we already know the 
entanglement of the corresponding pure state \ip(x,y)), 
we may accomplish our aim by the Vollbrecht- Werner 
technique Q, which gives the entanglement of p(x,y) in 
terms of that of \ip(x,y)} via the convex hull construc- 
tion: E p (x,y) = (co E^){x, y). Said in words, the en- 
tanglement surface z = E p (x,y) is the convex surface 
constructed from the surface z = E^,(x,y). 

The idea underlying the use of the convex hull is this. 
Due to its linearity in x and y, the state p{x, y) (|61|l can 
[except when (x, y) lies on the boundary] be decomposed 
into two parts: 



VSe^ 1 |GHZ) + Vy e i02 |W) + Vl-x-ye* |W). (65) 



p(x, y) =pp{x 1 ,y 1 ) + (1 - p)p(x2,V2), (69) 



. 2 



x 



i o 

FIG. 6: Entanglement of the pure state \%p(x, (1 — x)r)) = 

y/x |GHZ) + y/(l-x)r jW) + v'(l-x)(l-r)|W} vs. x and r. 
Note the symmetry of the surface with respect with r = 1/2. 
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FIG. 7: Entanglement of the pure states \ip(x,(l — x)r)} = 

y/x | GHZ) + y/(l -x)r | W) + v'(l-x)(l-r)|W} vs. r for var- 
ious values of x (from the bottom: 0.8, 0.85, 0.9, 0.92, 0.94, 
0.96, 0.98, 1). This reveals the nonconvexity in r for interme- 
diate values of x. 

with the weight p and end-points (xx,yx) and (£2,2/2) 
related by 

p xi + (1 - p).t 2 = x (70a) 
PS/i + (1 -p)V2 = y. (70b) 

Now, if it should happen that 

pE^xx.yx) + (1 -p)Ei jl (x2,y2) < E^(x,y) (71) 

then the entanglement, averaged over the end-points, 
would give a value lower than that at the interior point 
{x,y); this conforms with the convex-hull construction. 

It should be pointed out that the convex hull should 
be taken with respect to parameters on which the density 
matrix depends linearly, such as x and y in the example 
of p(x,y). Furthermore, in order to obtain the convex 
hull of a function, one needs to know the global struc- 
ture of the function — in the present case, E^(x,y). We 
note that numerical algorithms have been developed for 
constructing convex hulls [20J . 



FIG. 8: Entanglement of the pure states \ip\x, (1 — x ) r )) = 

Vi|GHZ} + y/(l -x)r|W) + V'(l-x)(l-r)|W) vs. x for 
various values of r (from the top: 0, 0.1, 0.2, 0.3, 0.5). This 
reveals the nonconvexity in x in the (approximate) interval 
[0.85,1]. 



As we have discussed, our route to establishing the en- 
tanglement of p(x,y) involves the analysis of the entan- 
glement of \tj}(x,y)), which we show in Fig.0] Although 
it is not obvious, the corresponding surface fails to be 
convex near the point (x,y) = (1,0), and; therefore in 
this region we must suitably convexify in order to obtain 
the entanglement of p(x,y). To illustrate the properties 
of the entanglement of \tp(x, y)) we show, in Fig. ^ the 
entanglement of \ijj(x,y)} along the line (x,y) = (0, s); 
evidently this is convex. By contrast, along the line 
x + 2y — 1 there is a region in which the entanglement is 
not convex, as Fig.[3]shows. The nonconvexity of the en- 
tanglement of \i/)(x, y)) complicates the calculation of the 
entanglement of p(x, y), as it necessitates a procedure for 
constructing the convex hull in the (as it happens, small) 
nonconvex region. Elsewhere in the xy plane the entan- 
glement of p(x, y) is given directly by the entanglement 
of \ip(x,y)). 

At worst, convexification would have to be undertaken 
numerically. However, in the present setting it turns out 
that one can determine the convex surface essentially an- 
alytically, by performing the necessary surgery on surface 
z = E^{x, y). To do this, we make use of the fact that if 
we parametrize y via (1 — x)r, i.e., we consider 

p(x,(l-x)r) = x|GHZ)(GHZ| + (l-x)r|W)(W| 

+(l-x)(l-r)\W)(W\, (72) 

where < r < 1 [and similarly for \?p(x,y))] then, as 
a function of (x, r), the entanglement will be symmetric 
with respect to r = 1/2, as Fig. makes evident. With 
this paramctrization, the nonconvex region of the entan- 
glement of |^) can more clearly be identified. 

To convexify this surface we adopt the following conve- 
nient strategy. First, we reparametrize the coordinates, 
exchanging y by (1 — x)r. Now, owing to the linearity, in 
r at fixed x and vice versa, of the coefficients x, (1 — x)r 
and (1 — x)(l — r) in Eq. (|72|) . it is certainly necessary 
for the entanglement of p to be a convex function of r at 
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FIG. 9: Entanglement of the mixed state p(x,y). 



fixed x and vice versa. Convexity is, however, not neces- 
sary in other directions in the [x, r) plane, owing to the 
nonlinearity of the the coefficients under simultaneous 
variations of x and r. Put more simply: convexity is not 
necessary throughout the [x, r) plane because straight 
lines in the (x, r) plane do not correspond to straight 
lines in the (x, y) plane (except along lines parallel either 
to the r or the x axis). Thus, our strategy will be to 
convexify in a restricted sense: first along lines parallel 
to the r axis and then along lines parallel to the x axis. 
Having done this, we shall check to see that no further 
convexification is necessary. 

For each x, we convexify the curve z = (x, (1 — x)r) 
as a function of r, and then generate a new surface by 
allowing x to vary. More specifically, the nonconvexity in 
this direction has the form of a symmetric pair of minima 
located on either side of a cusp, as shown in Fig. [7] Thus, 
to correct for it, we simply locate the minima and connect 
them by a straight line. 

What remains is to consider the issue of convexity 
along the x (i.e., at fixed r) direction for the surface just 
constructed. In this direction, nonconvexity occurs when 
x is, roughly speaking, greater than 0.8, as Fig. |S] sug- 
gests. In contrast with the case of nonconvexity at fixed 
r, this nonconvexity is due to an inflection point at which 
the second derivative vanishes. To correct for it, we lo- 
cate the point x = xq such that the tangent at x — xo is 
equal to that of the line between the point on the curve at 
xq and the end-point at x = 1, and connect them with a 
straight line. This furnishes us with a surface convexified 
with respect to x (at fixed r) and vice versa. 

Armed with this surface, we return to the (a;, y) 
parametrization, and ask whether or not it is fully con- 
vex (i.e., convex along straight lines connecting any pair 
of points). Said equivalently, we ask whether or not any 
further convexification is required. Although we have 
not proven it, on the basis of extensive numerical ex- 
ploration we are confident that the resulting surface is, 
indeed, convex. The resulting convex entanglement sur- 
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FIG. 10: Entanglement of the mixed state 

p{x,y=(l-x)/2) = x|GHZ)(GHZ| + l=s(|W)(W| + 

W)(W|) vs. x. Inset: enlargement of the region x £ [0.2, 0.3]. 
This contains the only point, (x,y) = (1/4,3/8), at which 
E p (x,y) vanishes. 

face for p(x,y) is shown in Fig. [3] Figure ITU1 exemplifies 
this convexity along the line 2y + x = 1. We have ob- 
served that for the case at hand it is adequate to correct 
for nonconvexity only in the x direction at fixed r. 



C. Comparison with the negativity 

This measure of entanglement is denned to be twice 
the absolute value of the sum of the negative eigenvalues 
of the partial transpose of the density matrix 0, ■ In 
the present setting, viz., the family p(x,y) of three-qubit 
states, the partial transpose may equivalently be taken 
with respect to any one of the three parties, owing to 
the invariance of p{x, y) under all permutations of the 
parties. Transposing with respect to the third party, one 
has 

N(p) =-2j2 A*, (73) 

Ai<0 

where the A's are the eigenvalues of the matrix p T3 , 

It is straightforward to calculate the negativity of 
p(x, y); the results are shown in Fig. 1111 Interestingly, for 
all allowed values of (x,y), the state p(x,y) has nonzero 
negativity, except at (x,y) = (1/4,3/8), at which the 
calculation of the GME shows that the density matrix 
is indeed separable. One also sees that /o(l/4,3/8) is a 
separable state from that fact that it can be obtained 
by applying the projection P4 to the (un-normalized) 

separable pure state (|0) + |1)) . The fact that the 
only positve-partial-transpose (PPT) state is separable 
is the statement that there are no entangled PPT states 
(i.e., no PPT bound entangled states) within this fam- 
ily of three-qubit mixed states. The negativity sur- 
face, Fig. 1111 is qualitatively — but not quantitatively — 
the same as that of GME. By inspecting the negativity 
and GME surfaces one can see that they present order- 
ing difficulties. This means that one can find pairs of 
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FIG. 11: Negativity of the mixed state p(x,y). 



states p(xi,y\) and p(x2,y2) that have respective nega- 
tivities Ni and N2 and GMEs E\ and E2 such that, say, 
Ni < N2 but Ei > £2. Said equivalently, the negativ- 
ity and the GME do not necessarily agree on which of 
a pair of states is the more entangled. For two qubit 
settings, such ordering difficulties do not show up for 
pure states but can for mixed states 0, ^2^. On the 
other hand, for three qubits, such ordering difficulties 
already show up for pure states, as the following exam- 
ple shows: A(GHZ) = 1 > N(W) = 2^2/3 whereas 
for the GME the order is reversed. We note, however, 
that for the relative entropy of entanglement En, one 
has E R {GUZ) = log 2 < E R {W) = log(9/4) 0, which 
for this particular case is in accord with the GME. 



VI. CONCLUDING REMARKS 

We have considered a rather general, geometrically mo- 
tivated, measure of entanglement, applicable to pure and 
mixed quantum states involving arbitrary numbers and 
structures of parties. In bipartite settings, this approach 
provides an alternative — and generally inequivalent — 
measure to the entanglement of formation. For multipar- 
tite settings, there is, to date, no explicit generalization 
of entanglement of formation |23|| . However, if such a 
generalization should emerge, and if it should be based 
on the convex hull construction (as it is in the bipartite 
case), then one may be able to calculate the entangle- 
ment of formation for the families of multipartite mixed 
states considered in the present Paper. 

As for explicit implementations, the geometric measure 
of entanglement yields analytic results in several bipartite 
cases for which the entanglement of formation is already 
known. These cases include: (i) arbitrary two-qubit 
mixed, (ii) generalized Werner, and (iii) isotropic states. 
Furthermore, we have obtained the geometric measure of 
entanglement for certain multipartite mixed states, such 
as mixtures of symmetric states. In addition, by mak- 



ing use of the geometric measure, we have addressed 
the entanglement of a rather general family of three- 
qubit mixed states analytically (up to root-finding) . This 
family consists of arbitrary mixtures of GHZ, W, and 
inverted-W states. To the best of our knowledge, cor- 
responding results have not, to date, been obtained for 
other measures of entanglement, such as entanglement 
of formation and relative entropy of entanglement. We 
have also obtained corresponding results for the nega- 
tivity measure of entanglement. Among other things, we 
have found that there are no PPT bound entangled states 
within this general family. 

A significant issue that we have not discussed is how 
to use the geometric measure to provide a classification 
of entanglement of various multipartite entangled states, 
even in the pure-state setting. For example, given a 
tripartite state, is all the entanglement associated with 
pairs of parts, or is some attributable only to the sys- 
tem as a whole? More generally, one can envisage all 
possible partitionings of the parties, and for each, com- 
pute the geometric measure of entanglement. This would 
provide a hierarchical characterization of the entangle- 
ment of states, more refined than the global characteri- 
zation discussed here. Another extension would involve 
augmenting the set of separable pure states with certain 
classes of entangled pure states, such as biseparable en- 
tangled, W-type and GHZ-type states |24|. 

Although there is no generally valid analytic procedure 
for computing the entanglement eigenvalue A max , one 
can give — and indeed we have given — analytical results 
for several elementary cases. Harder examples require 
computation, but often this is (by today's computational 
standards) trivial. We note that in order to find A max 
for the state it is not necessary to solve the nonlinear 
eigenproblem JSJ ; one can instead appropriately parame- 
terize the family of separable states \4>) and then directly 
maximize their overlap with the entangled state i.e., 
A max = max0 1 1 {(75|-0) 1 1 ■ As an aside, we mention there 
exist numerical techniques for determing E-p (see, e.g., 
Ref. jl^l). We believe that numerical techniques for solv- 
ing the geometric measure of entanglement for general 
multipartite mixed states can readily be developed. 

The motivation for constructing the measure discussed 
in the present Paper is that we wish to address the degree 
of entanglement from a geometric viewpoint, regardless 
of the number of parties. Although the construction is 
purely geometric, we have related this measure to entan- 
glement witnesses, which can in principle be measured 
locally ^^|. Moreover, the geometric measure of entan- 
glement is related to the probability of preparing a sin- 
gle copy of a two-qubit mixed state from a certain pure 
state |l8j. Yet it is still desirable to see whether, in gen- 
eral, this measure can be associated with any physical 
process in quantum information, as are the entanglement 
of formation and distillation. 

There are further issues that remain to be explored, 
such as additivity and ordering. The present form of 
entanglement for pure states, E s - m 2 = 1 — A 2 , is not ad- 
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ditive. However, one can consider a related form, E\ og = 
— log A 2 , which, e.g., is additive for \iP}ab ® \*P)cd, i-e., 

Elog(\i>l)AB® \^2)cd) = ^log ( | "01 ) Afi) + ^log {\^2)cd) ■ 

(74) 

This suggests that it is more appropriate to use this loga- 
rithmic form of entanglement to discuss additivity issues. 
However, it remains to check whether it is an entangle- 
ment monotone when extended to mixed states by convex 
hull. 

As regards the ordering issue, we first mention a result 
of bipartite entanglement measures, due to Virmani and 
Plenio .2j| , which states that any two measures with con- 
tinuity that give the same value as the entanglement of 
formation for pure states are either "identical or induce 
different orderings in general'.' This result points out that 
different entanglement measures will inevitably induce 
different orderings if they are inequivalent. This result 
might still hold for multipartite settings, despite their 
discussion being based on the existence of entanglement 
of formation and distillation, which have not been gener- 
alized to multipartite settings. Although the geometric 
measure gives the same ordering as the entanglement of 
formation for two-qubit mixed states [see Eq. H45|) ]. we 
belive that the geometric measure will, in general, give 
a different ordering. However, it is not our intention to 
discuss the ordering difficulty in the present Paper. Nev- 
ertheless, it is interesting to point out that for bipartite 
systems, even though the relative entropy of entangle- 
ment coincides with entanglement of formation for pure 
states, they can give different orderings for mixed states, 
as pointed out by Verstraete et al. |22| . 

We conclude by remarking that the measure discussed 
in the present Paper is not included among the infinitely 
many different measures proposed by Vedral et al. [26j. 
Those measures are based on the minimal distance be- 
tween the entangled mixed state and the set of separable 
mixed states. By contrast, the measure considered here is 
based upon the minimal distance between the entangled 
pure state and the set of separable pure states, and it is 
extended to mixed states by a convex hull construction. 
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APPENDIX A: THE VOLLBRECHT- WERNER 
TECHNIQUE 

In this Appendix, we now briefly review a technique 
developed by Vollbrecht and Werner [(| for computing 
the entanglement of formation for the generalized Werner 
states; this turns out to be applicable to the computation 
of the sought quantity E sin 2 . We start by fixing some 
notation. Let 

(a) K be a compact convex set (e.g., a set of states 
that includes both pure and mixed ones); 

(b) M be a convex subset of K (e.g., set of pure states); 

(c) E : M — ► i? U {+00} be a function that maps ele- 
ments of M to the real numbers (e.g., E = E sin 2); 
and 

(d) G be a compact group of symmetries, acting on K 
(e.g., the group U <g> U^) as a g : K — > K (where 
a g is the representation of the element g £ G) that 
preserve convex combinations. 

We assume that a g M C M (e.g., pures states are 
mapped into pure states), and that E(a g m) = E(m) 
for all to £ M and g £ G (e.g., that the entanglement of 
a pure state is preserved under a g ). We denote by P the 
invariant projection operator defined via 

Pk = J dga g (k), (Al) 

where k £ K. Examples of P are the operations Pi and 
P2 in the main text. Vollbrecht and Werner also defined 
the following real-valued function e on the invariant sub- 
set PK: 

e(x) = inf {E(m)\m £ M, Pm = x} . (A2) 

They then showed that, for x £ PK, 

coE(x) = coe(x), (A3) 

and provided the following recipe for computing the func- 
tion co E for G-invariant states: 

1. For every invariant state p (i.e., obeying p = Pp), 
find the set M p of pure states a such that Pa = p. 

2. Compute e(p) = inf {E(a)\a £ M p }. 

3. Then coE is the convex hull of this function e. 
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